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A  CONTINUUM  DIFFUSION  MODEL  FOR  VISCOELASTIC  MATERIALS 


By 

Y.  Weitsman 

Mechanics  and  Materials  Center,  Civil  Engineering  Department 
Texas  A&M  University,  College  Station,  Texas  77843-3136 

Abstract 

A  model  for  diffusion  in  polymers  is  established  from  basic  principles  of 
irreversible  thermodynamics,  employing  the  methodology  of  continuum  me¬ 
chanics.  The  polymeric  materials  are  considered  to  respond  viscoelastically, 
with  ageing. 

It  is  shown  that  effects  of  stress  on  diffusion  and  certain  anomalies  in 
the  moisture  sorption  process  can  be  explained  by  the  present  model. 
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1.  Introduction 

The  processes  of  moisture  transport  and  sorption  in  polymers  have  been 
studied  by  numerous  investigators  for  more  than  fifty  years.  Most  of  these 
studies,  which  involved  experimental,  analytical,  and  materials  science 
aspects,  were  conducted  by  researchers  in  the  fields  of  physical  and  polymer 
chemistry.  These  investigations  were  mostly  focused  towards  applications  in 
membrane  technology  and  therefore  concerned  aspects  of  permeability  and 
seepage.  It  is  far  beyond  the  scope  of  the  present  article  to  list  the  many 
significant  contributions  made  to  this  field.  References  !1)  -  [9]  list  some 
of  the  many  outstanding  review  articles  on  the  subject. 

During  the  past  fifteen  years  the  transport  of  moisture  in  polymers 
raised  concern  and  interest  in  composite  materials  technology.  Moisture  plays 
a  significant  role  in  affecting  the  mechanical  response  of  fiber-reinforced 
polymeric  composites  and  in  influencing  the  behavior  of  adhesive  joints  which 
are  increasingly  employed  in  composite  structures.  These  concerns,  which 
focus  mainly  on  mechanical  effects,  produced  more  than  500  publications  on  the 
subject.  A  review  and  summary  of  some  of  these  works  is  given  in  references 
[10]  and  [11|. 

Inspite  of  many  overlaps  the  investigations  in  the  two  foregoing 
categories  reveal  distinctions  that  reflect  the  disparate  scientific 
disciplines  from  which  they  derive  -  polymer  science  in  one  case  and  applied 
mechanics  in  the  other.  It  may  thus  be  worthwhile  to  present  an  approach  to 
theory  of  diffusion  that  employs  the  methodology  of  continuum  mechanics  and 
its  well  developed  constitutive  formalism.  This  approach  is  especially  useful 
in  the  modelling  of  coupling  effects,  such  as  between  stress  and  diffusion,  as 
detailed  in  a  recent  work^^^L 

Tf  should  be  noted  that  the  most  crucial  issue  associated  with  diffusing 
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substances  in  polymeric  composites  is  their  detrimental  effects  on  structural 
strength,  integrity  and  durability.  In  certain  circumstances,  diffusing 
substances  were  noted  to  damage  the  composites  by  debonding  the  fiber/matrix 
interfaces.  This  issue  is  a  topic  of  current  research  and  its  modelling  is 
hampered  by  uncertainties  regarding  the  damage  mechanisms.  Several 
investigations  concerning  moisture  induced  damage  in  composites  were  reviewed 
in  ref.  [11]  and  an  attempt  at  a  continuum-level  modelling  of  the  phenomenon 
is  provided  in  ref.  [13] .  The  issue  of  damage  will  not  be  discussed  in  the 
present  paper. 

2.  Basic  Equations 

Consider  a  solid  body  B  occupying  a  material  volume  V  bounded  by  a 
surface  A.  Let  the  solid,  of  mass  density  p^,  absorb  vapor  through  its 
boundary  and  let  m  denote  the  vapor-mass  per  unit  volume  of  the  solid.  Also, 
let  X  be  the  position  of  a  solid  mass  particle  in  the  deformed  configuration 
that  corresponds  to  the  place  X  in  the  undeformed  state,  and  let  f,  g  and  v 
denote  fluxes  of  vapor -mass  and  of  heat,  and  the  velocity  of  the  solid 
particles,  respectively. 

In  addition,  let  u  and  s  be  the  internal  energy  and  entropy  densities  of 
the  solid/vapor  mixture  per  unit  solid  mass,  and  let  o.j  and  T  denote  the 
components  of  the  Cauchy  stress  due  to  mechanically  applied  loads,  and 
temperature,  respectively. 

A  proper  accounting  of  the  state  of  the  solid/vapor  mixture,  which  is  a 
thermodynamically  open  system,  is  obtained  by  considering  each  element  in 
thermodynamic  equilibrium  with  a  reservoir  containing  vapor  at  pressure  p, 
density  p,  and  internal  energy  and  entropy  densities  u  and  s  respectively 
114],  115],  116). 

Conservation  of  the  solid  and  vapor  masses  gives 
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*s  ^  "s  ^  ^  °  (1) 

m  =  -V  •  f  (2) 


Conservation  of  energy  over  B  reads 

ST  i  "s  "  ‘  i  ■  I 

-  1  p  —  n.dA  -  J  uf.n.  dA  (3) 

A  0  A  '  ' 

The  third  integral  on  the  right  side  of  (3)  expresses  the  mechanical 
power  due  to  vapor  flux,  observing  that  f^/o  corresponds  to  vapor  velocity. 
The  last  integral  in  (3)  expresses  the  rate  of  vapor-borne  energy. 

The  entropy  inequality  reads 

^  ;  0  S  dV  >  ;  -(q./T)n.  dA  -  I  sf^n  dA  (4) 

Where  the  last  integral  in  (4)  expresses  the  rate  of  vapor-borne  entropy. 
Application  of  Green’s  theorem  to  (3)  and  (4),  and  employment  of  (2), 
yields^ 


and 


(5) 


^s^  ^  -  ■‘^i.i  +  T  s  m  (6) 

where  h  =  (p/o)  +  u  is  the  enthalpy  of  the  vapor  in  the  hypothetical  reservoir 
and  =  T^^. 

Elimination  of  q^j  between  (5)  and  (6)  yields  the  following  expression 
for  the  "reduced  entropy  inequal ity" ^ 
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-p 


y  ^so  ij 


so"  "SO'  'ij  'IJ  '■’V  '  =  1  ~  -  s  i 


(7) 


In  (7)  ?  =  u-Ts  -  7 —  'ii^ii  Gibbs  free  energy  and  p  =  h-Ts  is 
"so 

the  chemical  potential  of  the  vapor  in  the  hypothetical  reservoir.  In  eqn. 
(7)  0^^  denotes  the  initial  mass  density  of  the  solid  and  infinitesimal 
deformation  has  been  assumed,  whereby  «  1. 


3.  Ageing  Viscoelastic  Response  With  Moisture. 

The  viscoelastic  behavior  of  high  polymers  has  been  associated  with  the 
various  degrees  of  freedom  afforded  by  the  motions  of  the  long  chain  molecules 
!17),  1181,  |19)  normal  modes  of  these  motions  can  be  represented  by  N 
scalar-valued  internal  state  variables  (r  =  1,  ...  N).  This 
representation  motivated  subsequent  researchers^ to  derive  a  theory  of 
viscoelasticity  from  fundamental  concepts  of  irreversible  thermodynamics,  with 
irreversibility  stemming  from  the  foregoing  variables  through  "growth 
laws"  that  correspond  to  viscous  resistance. 

The  physical  ageing  of  polymers  is  generally  ascribed^  ^  to  the 
introduction  of  a  "free  volume"  into  those  materials  upon  their  cooling  across 
the  glass  transition  temperature  Tg.  The  free  volume  is  a  non-equilibrium 
thermodynamic  internal  state  variable,  which  decays  spontaneously  with  the 
passage  of  time  and  influences  the  rate  of  all  time-dependent  processes  within 
the  polymer. 

Consider  therefore 


$  =  (ji(p^ j ,  m,  T,  Y^,  v^)  (8) 

where  v^  denotes  the  non-dimensional  free-volume  fraction. 

Since  eqn.  (7)  cannot  be  violated  by  any  process,  the  methodology  of 
continuum  mechanics^ ^  yields  the  following  constitutive  equations 
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s  = 


-  -D 


SO  3C.  . 


3  f 

3T 


(9) 


and 


- 

”  “^so  sm 


-R^  -r^  -  R  -  (q^/T)g.  -  S  g.f.  >  0 


(10) 


In  eqn.  (10),  the  affinities  and  R  are  defined  by 


=  0, 


r  f 


The  aging  process  is  associated  with  the  spontaneous  collapse  of  the  free 
volume  1251,127)  initial  value  v^q,  which  is  created  by  the  cooling 

of  a  polymer  across  is  glass  transition  temperature  Tg  down  to  a  temperature 
Tq.  In  addition,  the  subjection  of  a  polymer  externally  to  m,  or 

aT  =  ^'^0*  trigger  an  irreversible  process  which  will  cause  the  internal 

p 

variables  to  drift  spontaneously  toward  their  equilibrium  values  •  The 
two  abovementioned  processes  are  assumed  independent  of  each  other 

1251,127]^  Consequently,  while  both  and  v^  in  eqn.  (10)  denote  derivatives 

with  respect  to  time,  the  time-lapse  t  for  is  not  necessarily  equal  to  the 

time-span  t  for  v^. 

Assume  'r^,  v^  «  1.  A  Taylor  expansion  about  the  initial  state 

(y^  =  0,  v^  =  v^q)  then  reads 


^o"  =  ^so 


(Jj  +  $  Y  +  $  (y^-  V-  )  +  i  If  _Y..Y.,  o’  ♦  (Vf-  Vr„) 

0  r  r  v'  f  fo'  2  rq  r  q  2  vv'  f  fo' 


+  H .  0 .  T . 


(12) 


In  eqn.  (12)  =  <fQ(o^j,  fn.T,  0,  and  i)>^,  4.^^,  4^^  denote  partial 
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derivatives  at  =  0,  v^  =  v^^,  all  of  which  depend  on  o-j,  m,  and  T. 

According  to  the  extremum  principles  of  thermodynamics  4>  attains  its 
minimum  at  the  equilibrium  values  and  v^  =  Vp  Therefore,  at 
equilibrium  we  have 


1^-0,  ^~0.  and 


3y 


3V 


,2  ,2  , 

2  17^7“  ^’^r^'^s  2  - r  ^  ° 

Expanding  4)  about  y^  and  v®  =  0  we  have 

^0^  "  ^o'^e  1  Vq  (V  (^q-  ^q)  1  Sv  ''f  >  (1®) 


In  view  of  the  assumption  that  y^.  y^.  v^,  v^^  «  1,  expansions  (12)  and 
(13)  are  equal,  implying  in  particular  that 


♦  vv  ”  '^vv  *  ^rq  ~  '*’rq 


Therefore,  >  0  and  iti^^  are  components  of  a  symmetric,  positive 
definite  matrix. 

Employing  the  familiar  assumption  of  viscous-like  resistance^ let 

dv. 


and 


R  =  -b(m,  T,  Vf) 


=  -b  (m,  T,  V,) 
r  rq'  ’  f' 


dt 


(14) 

(15) 


By  Onsager's  principle  b^q  =  bq^  and,  in  view  of  the  dissipation 
inequality  (10),  b>0  and  b^q  are  components  of  a  semi  positive  definite, 
symmetric  matrix.  As  noted  earlier,  the  time  t  and  t  in  eqns.  (14)  and  (15) 
are  distinct  because  the  quench-time  t  usually  occurs  earlier  than  the 
exposure-time  t,  namely  t  =  T-tg. 

Equations  (11),  (13)  and  (14)  give 
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whereby 


e 

p .  $ 
sovv 


with  d;  =  d-. /b 
Equations  (11), 


and 


^-1  ^  ^  e 

■  V  "'so'^vv 


(12)  and  (15)  yield 


(16) 

(17) 


0  $  y  +  b  y  ““0  c 

SO  rq’^q  rq  ^q  r 


(18) 


Since  and  b^q  are  symmetric,  positive  definite  matrices*  they  can  be 
diagonalized  simultaneously  and  eqn.  (18)  is  expressible  in  the  form 

(p  =  1,...N.  no  sum  on  p)  (19) 


A  n  ^ 

$  >  +  B  Y 

SO  p  'p  p  p 


-  p  $ 

so*p 


Note  that  bpq,  and  therefore  also  Bp,  depend  on  the  quency-time  i  as  well 
as  on  the  exposure-time  t  through  their  dependence  on  m,  T  and  in  eqn. 
(15). 

Following  Schapery  1221-1241  355^1^,0  t:hat  all  :t>^q  have  a  common  dependence 
on  o^j,  m,  and  T  namely  =  ^G^^i  j’”’’^^‘’’rq  constant  and,  simi  larly, 
b^q  =  a(m,T,Vf)  0^°.  These  common  "shift-factors"  aQ  and  a  carry  over  to 
eqn.  (19)  whose  solution,  subject  to  the  conditions  y^(0)  =  0,  Yj.(“)  = 
reads 


A 


=  -K 


.(1-e  """'P) 


(20) 


♦The  case  of  a  semi-positive  definite  b^g  is  omitted  here  for  the  sake  of 
brevity.  This  case  corresponds  to  a  viscoelastic  fluid,  as  modelled  by  a 
free  dash  pot  attached  to  Kelvin  or  Maxwell  elements. 
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where 


<  =  - ^ 

p 


= 


and 


G"P 


P  0  <t 
^  SO  p 


(no  sum  on  p,  p  =  1,2,. ..N) 


(21) 

(22) 


Aiming  at  linear  stress-strain  behavior,  consider  o^j/E  «  1  where  E 
denotes  a  characteristic  modulus.  To  retain  linearity  assume  that  all 
retardation  times  do  not  depend  on  o.j.  Hence,  in  eqns.  (17)  and  (20), 
\  =  T^(m,T,v^)  and  =  Tp(m,T,v^)  p  =  1,...N.  Furthermore,  to  establish 
the  linear  stress-strain  relations,  consider  the  expansion  in  eqn.  (12) 
truncated  after  three-terms,  which,  resorting  to  the  diagonalized  form,  reads 


C  $  =  0  1$  +  itY  (v,  -  V-_)  I 

so  so  0  r'r  ^v'  f  fo' 


(23) 


Eqn.  (9)|  gives 


^ij  ‘  ■‘^SoUo.j  3o^j  V  3o^j  ^''f  ■  ''fo^  * 


(24) 


The  specific  forms  which  follow  from  eqn.  (24)  depend  on  the  material 
symmetry  at  hand.  Considering  isotropy  we  have 

0(a.y  m,  T,  Y^,  v^)  =  ^ij'^ij’  ^r’  ''f^ 

which,  for  linear  behavior  gives  the  well  known  forms 

“so*o  "  -  "o("-''>“k.»,k 

“so*r  *  ■*  '-r<"-''>”kk  *  "r<”-’>°kk»..  * 

It  is  advantageous  to  express  Lq,  and  in  terms  of  expansional 
coefficients  a  and  b,  namely 


♦The  stress  invariant  is  inadmissible  for  linear  behavior. 
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It  can  be  noted  from  eqns.  (20)  and  (21)  that  Yp  depends  on  through 
the  presence  of  $p  in  Kp.  However,  for  sufficiently  small  stresses  Kp 
themselves  can  be  expanded  according  to  (25),  whereby  upon  collecting  like- 
powers  in  a.j,  the  form  of  the  strain-stress  relationship  given  in  eqn.  (27) 

A 

remains  valid  with  stress- independent  y^. 

The  stress  strain  relations  (27)  are  of  the  familiar  form  employed  in 
linear  viscoelasticity,  except  that  the  retardation  spectra  incorporate  now 


ageing  effects  and  all  instantaneous  compliances  age  with  time.  Recalling  eqn 
(20),  for  a  discrete  spectrum  of  retardation-times  expressions  (27)  read 

.  =  (a^  ^  a,(l-e  '^)^Ia^(l-e  '')  IcT 

-t-  u  +  6,(l-e  )  +  L  s  (1-e  )lm 

+  1B  +B,(l-e  '')  +  r  B  (1-e  )lo..  6.. 

°  ^  -c/:  r  -c/t^ 

+  |J  +  J^(l-e  ^  I  J  (1-e  ’^)lo..  (28) 


In  eqn.  (28),  the  quantities  a,  e,  B  and  J  depend  on  m  and  T,  the 
"reduced  times"  t  and  c  both  depend  on  m,  T  and  the  ageing  time  t.  However, 
in  contrast  to  ;,  the  reduced  time  t  ’S  measured  from  the  time  of  application 
of  o.j,  m,  and  T.  For  fluctuating  environments  and  stresses,  expression  (28) 
takes  the  form  of  convolution  integrals.  The  case  of  a  continuous  spectrum  of 
retardation  times  can  be  represented  in  the  familiar  manner  by  means  of 
retardation  integrals. 
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Stress  Assisted  Diffusion  in  Ageing  Viscoelastic  Media 


Employing  eqns.  (9)3,  (23),  (25)  and  (26)  we  obtain 

3A^  3A^(ni,T)  3A^(ni,T) 

^  ~  3in  3m  V  3m  ^''f‘  '^fo^ 


3a_  3a 


^  liT  ^  V  ^  iiT  (V-  ''fo^’^'^ 

3(mB^)  3(nie^) 


3(mB^)  3(nie^) 

^  3m  sin  ^^f'  '^fo^^'^kk 


3M^  3M^  3M^ 

^  I-  iiT  ^  iF  V  ^  IJT  (''f-  "fo^^^kk  °s 


oNq  3N^  3N^ 

^  ^ '  3m  ^  3m  ''r  ^  3m  ^'^f~  '^fo^^^ 


k4°!ik 


The  boundary  conditions  for  moisture  transport  are 

-(x,t)  =  ^^(x,t)  X  on  boundary  (30) 

where  represents  the  chemical  potential  of  the  ambient  vapor. 

Assume  for  simplicity  that  the  expansional  coefficients  a  and  a  do  not  vary 
with  m.  In  this  case  eqns.  (29)  and  (30),  together  with  eqn.  (20),  yield 

3A  -^/*r  -  t/i 

-^  =  +  z  W^(m,T)  (1-e  *  W(m.T)  (1-e  ") 

r 

-t/'  -c/i 

+  1-6  +  z  6^(l-e  )  +  6y(l-e  )l°kk 

r 

+  l-Y(m,T)  +  I  Y^(m,T)  (1-e  '")  +  Y^(l-e  '^)l°kk°aii 

r 

+  l-Z(m.T)  +  I  Z^(m.T)  (1-e  ^^‘'■)  +  Z^(l-e  1  Okj,Ok£ 

Eqn.  (31)  illuminates  the  effects  of  viscoelastic  retardation  and  ageing, 
as  well  as  of  stress,  on  the  diffusion  process.  In  the  absence  of  those 

effects  the  boundary  condition,  which  reads  3A^/3m  =  translates  into  the 
familiar  statement  m(x,t)  =  mg(t),  x  on  boundary  and  m^  prescribed.  In 
contrast,  eqn.  (31)  states  that  the  equilibrium  boundary  value  is  approached 
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gradually  with  time  (even  for  exposure  to  constant  ambient  vapor  pressure),  it 
is  affected  by  age,  and  depends  quadratically  on  the  applied  stress. 

Consider  now  the  diffusion  coefficient  D.  Since  this  coefficient 
reflects  micro-level  phenomena,  its  form  cannot  be  derived  from  a  continuum 
model*.  Molecular-level  considerations  of  diffusion  in  solids^^^^  and 
fluids^^^^  suggest  that 

0  -  Vacancy  in  host  material  •  Agitation  energy  of  diffusing  substance 
In  the  present  case  it  is  reasonable  to  expect  that  the  vacancy  would 
depend  mostly  on  v^  and  secondarily  on  m,T  and  (through  their  influence  on 
swelling  strains).  On  the  other  hand,  the  agitation  energy  would  depend 
mostly  on  T. 

Accordingly 


0  -  0Q(m,T,v^(c(T)),  J*  F(t(t)  -  Ur))  d-l  .  (exp(-E/RT))  (32) 

0 

In  the  absence  of  a  satisfactory  molecular  theory  for  the  glassy  state  at 
the  present  time  it  appears  expedient  to  base  the  form  of  0  on  empirical 
evidence. 


In  isotropic  materials,  the  flux  of  moisture  at  fixed  and  T  is  taken 


to  be** 


f.  =  -D  — 

1  3X. 


Consequently,  eqn.  (2)  yields  the  field  equation 


(33) 


*The  principle  of  equi-presence  suggests  that  D  =  D(o^  . ,m,T,Y^, v^)  but 
provides  no  further  insight.  ^ 

**When  m,T  and  vary  in  space  couplings  occur  between  the  fluxes  of 
moisture  and  of  temperature  and  both  depend  on  the  gradients  of  p 
and  of  the  invariants  of  stress 


In  view  of  eqns.  (29)  and  (32),  relations  (33)  and  (34)  are  extremely 
complicated  even  at  T  =  Tq  and  the  complexity  is  due  to  the 

presence  of  m  in  D  and  v^  (through  ^(t),  and  c(t)  in  the  compliances  listed  in 
eq.  (29)).  However,  inspite  of  the  cumbersome  details  and  the  paucity  in 
information  several  features  of  the  form  of  equation  (34)  emerge. 
Accordingly,  the  diffusion  equation: 

(1)  Contains  non-linear  terms  in  the  moisture  gradient  am/sx^. 

(2)  Follows  a  time-retardation  process  akin  to  mechanical  viscoelastic  response. 

(3)  Varies  non-linearly  with  external  stresses 

(4)  Exhibits  an  ageing  behavior  characteristic  of  glassy  polymers. 

5.  Two-Phase  Diffusion 

The  non-Fickean  weight  uptake  of  water  vapor  by  epoxy  resins  can  be 
described  by  a  two-phase  diffusion  model.  This  model  assumes  that  the  vapor 
within  the  polymer  divides  into  two  portions  -  a  bound  phase  m^.  and  a  mobile 
phase  m^.  Details  are  given  in  refs.  [331,  134],  where  the  resin's 

viscoelasticity  has  been  discarded. 

It  is  interesting  to  note  that  the  above  model  can  be  derived  in  the 
context  of  the  thermodynamic  theory  developed  in  this  work.  In  this  case,  let 

=  ^0^  ('’ij*  "’a*  "’c*  ""r’  ''f) 

In  analogy  with  sections  2-4  we  now  have 

a*  3* 

'^a  so  am^  ’  ^c  so  am^ 

d  C 

Equation  (29)  and  the  boundary  condition  (31)  pertain  now  to  the  chemical 
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potential  of  the  mobile  phase  u,,  so  that  all  partial  derivatives  in  those 

a 

equations  should  be  taken  with  respect  to  riig^  instead  of  m. 

Similarly,  the  flux-gradient  relation  now  involves  the  mobile  phase  only, 
so  that  eqn.  (33)  is  replaced  by 

3  u  ^ 


f.  =  -D  — ^ 

1  3X^ 


(36) 


Since  the  balance  equation  (2)  remains  valid,  equation  (34)  is  replaced 
now  by 


a  (m  +  m  )  .  3^ 

- =  -£-  (0  — ^) 

3t  3X^  '  3X.^ 


(37) 


The  system  of  equations  is  supplemented  by  the  phase  interaction 

relationl33l.l34| 


3m 

TT  '  ’’"c  -  ‘’"a 


(38) 


where  P  and  Q  are  empirical  factors. 

The  dependence  of  field  equations  and  boundary  conditions  on  retardation- 
times,  physical  ageing  and  stress  is  retained  for  the  case  of  two-phase 
diffusion  and  is  expressible  in  the  same  forms  as  presented  in  section  4. 


6.  Experimental  Observations. 

The  present  constitutive  model  can  be  employed  to  explain  certain 
experimental  observations  on  diffusion  in  polymers,  as  shown  in  Figs.  1-9 
below. 

Fig.  1  exhibits  the  uptake  of  water  vapor,  from  an  ambient  environment 
maintained  at  a  fixed  level  of  97%  R.H.,  by  unidirectionally  reinforced 
graphite/epoxy  coupons.  The  per-cent  weight  gain  is  plotted  vs.  v/time.  Note 
the  discrepancy  between  data  and  predictions  based  upon  the  linear  Fickean 


model . 


On  the  other  hand,  when  the  retardation  process  the  epoxy  is  taken  in 
consideration  by  means  of  a  boundary  condition  of  the  form 

-tA, 

m  (boundary)  =  m,  +  m,  (1-e  )  (39) 

it  is  possible  to  obtain  an  excellent  fit  between  data  and  theory.  In  eqn. 
(39)  m^  corresponds  to  the  instantaneous  response  and  mj^  is  the  contribution 
of  polymeric  retardation  processes.  The  corresponding  prediction,  with  m^  = 
1.2%,  m^  =  0.467%,  and  =  2.63x10^  min.,  is  drawn  with  a  solid  line  in  Fig. 
1.*  It  is  worth  noting  that  the  above  value  of  falls  roughly  in  the  middle 
of  the  spectrum  of  mechanical  retardation-times  for  epoxy. 

Another  form  of  history-dependence  of  moisture  sorption  in  epoxy  is 
suggested  in  Fig.  2^^^^  According  to  those  data  there  exists  a  distinct 
difference  in  the  moisture  up-take  of  "as  prepared"  samples  and  the  re¬ 
absorption  by  saturated-and-dried  samples.  This  behavior  can  be  attributed  to 
the  possibility  that  re-absorption  may  have  occurred  before  the  complete 
recovery  of  the  retardation  processes  caused  by  the  earlier  exposure. 

It  is  possible  to  match  data  of  the  type  shown  in  Figure  2  with  the 
boundary  condition  expressed  in  eqn.  (39).  The  incomplete  viscoelastic 
recovery  at  the  onset  of  the  second  exposure  to  moisture  can  be  introduced  by 
an  enhanced  value  of  m^,  compensated  by  a  reduced  value  of  m^,  maintaining  m^^ 
+  mK  =  m  =  constant.  Alternately,  the  incomplete  viscoelastic  recovery  may 
be  expressed  by  an  increased  value  of  t^.  The  trends  are  exhibited  in  Figs.  3 
and  4  respectively. 

Effects  of  stress  on  diffusion  are  shown  in  Figs.  Fig.  5 

shows  the  per-cent  weight  gain  of  water  vapor,  maintained  at  97%  R.H.,  by 

♦An  analytic  solution  for  boundary  condition  (39)  is  available  in  Ref.  (35). 
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epoxy  coupons  that  were  subjected  to  various  levels  of  uniaxial  tensile 
stress.  The  weight-gain  is  plotted  vs. ^ time  and  data  are  compared  against 
predictions  of  the  linear  Fickean  law,  showing  discrepancies  that  increase 
with  stress.  The  maximal  level  of  moisture  up-take  reached  during  the  tests 
(which  lasted  about  seven  weeks)  are  plotted  vs.  stress  in  Fig.  6.  Although 
the  above  maximal  levels  of  moisture  uptake  are  not  the  equilibrium  saturation 
levels.  Fig.  6  shows  a  trend  which  suggests  a  quadratic  relation  between 
saturation  levels  and  stress,  as  suggested  in  section  4. 

In  Figs.  5  and  6  the  ultimate  tensile  stress  was  o^-j^  =  5.17  MPa. 

The  dependence  of  water  absorption  by  polymers  on  compressive  stress 
levels  was  reported  in  ref.  (38|.  In  particular,  the  sorption  of  distilled 
water  by  epoxy  coupons  subjected  to  uniaxial  compression  is  shown  in  Fig.  7. 
Inspite  of  the  relatively  short  duration  of  these  tests,  which  were  terminated 
much  below  saturation,  they  indicate  a  non-linear  relationship  between  stress 
and  saturation  levels. 

A  similar,  though  less  direct,  non-linear  relation  between  stress  and 
saturation  levels  is  inferred  from  Fig.  This  figure  shows  the  reduced 
permeability  associated  with  the  flow  of  oxygen  gas  across  a  membrane  made  of 
extruded  PET,  when  the  membrane  is  subjected  to  various  levels  of  uniaxial  and 
biaxial  tensile  stresses. 

A  more  intricate  kind  of  sorption  behavior  is  shown  in  Fig.  gl^^i.  Since 
earlier  experiments  have  shown  that  equilibrium  solubility  decrease  with 
temperature^ ,  one  might  expect  to  observe  water  desorption  when  a  specimen 
equilibrated  at  75°C  is  transferred  to  a  water  bath  at  However,  Fig.  9 
shows  that  the  opposite  takes  place.  A  similar  phenomenon  was  also  observed 
in  graphite/epoxy  compos ites^^^L 

One  possible  explanation  of  this  behavior,  in  addition  to  those  suggested 


in  [361,  is  that  the  slow,  decelerating  collapse  of  the  free  volume  (initiated 
when  the  sample  was  c/eated)  is  reactivated  in  the  saturated  epoxy  when  it  is 
cooled  from  /S'-C  down  to  23''C.  Indeed,  it  was  indicated  in  1361,  that  for 


saturated  epoxy  Tg  is  lowered  to  lOS^C.  According  to  eqns.  (29)  and  (30)  the 
rejuvination  of  the  process  of  free  volume  collapse  would  indeed  be  reflected 
in  an  increased  saturation  level. 

6.  Concluding  Remarks 

In  this  article  it  was  shown  that  the  constitutive  formalism  of  continuum 
mechanics  can  be  employed  to  model  several  aspects  of  the  process  of  diffusion 
in  polymers.  Certain  micro-level  characteristics  of  the  polymer  structure, 
such  as  free  volume  and  motion  degrees-of-freedom,  were  incorporated  as 
internal  state  variables.  The  constitutive  formalism  led  to  a  linear 
viscoelastic  model  with  ageing  and  resulted  in  diffusion  relations  that  were 
coupled  with  stress  and  with  the  time-dependent  material  response. 

On  the  other  hand,  it  was  shown  that  the  continuum  mechanics  formalism 
provides  only  scant  information  about  material  parameters  such  as  the 
diffusion  coefficient.  This  shortcoming  is  attributable  to  the  fact  that 
continuum  mechanics,  which  aims  at  the  modelling  of  macro-level  phenomena,  is 
unsuitable  for  describing  inherently  micro-level  aspects.  This  conclusion 
indicates  that  an  interdisciplinary  effort,  combining  continuum  mechanics  and 
polymer  science,  is  required  for  a  comprehensive  formulation  of  the  diffusion 
process. 
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Figure  Titles 

Fig.  1:  Moisture  weight  gain  (%)  in  a  12-ply  unidirectional  AS4/3502 
graphite/epoxy  laminate  under  exposure  to  constant  ambient  relative 


humidity  of  97?^,  Vs.  /time.  Temperature  =  73°C.  Data  (dots), 
predictions  of  Fick's  law  (dashed  line),  and  of  eqn.  (39)  (solid 
line).  Values  of  m^,  m^^  and  given  in  text. 

Fig.  2:  Solubility  S  (in  %)  Vs.  non-dimensional  /time  of  epoxy  samples 
immersed  in  water  at  23“C.  Open  circles  -  samples  as  prepared,  solid 
circles  -  samples  previously  swollen  and  dried.  After  Apicella  et. 
al . ,  reference  1361. 

Fig.  3:  Moisture  weight  gain  (%)  Vs.\/time  for  test  conditions  as  in  Fig.  1 

and  predictions  according  to  eqn.  (39).  m^  and  mj^  as  in  text  (dashed 
line)  and  m^  =  1.25  m^,  -  0.25  m^  (solid  line).  In  both 

cases  ^  is  same  as  in  text. 


Fig.  4:  Moisture  weight  gain  (% )  Vs.  i/time  for  test  conditions  as  in  Fig.  1 
and  predictions  according  to  eqn.  (39).  as  in  text  (solid  line) 
and  =  t^/2  (dashed  line).  In  both  cases  m^  and  mj^  are  same  as  in 
text. 

Fig.  5:  Average  values  of  moisture  weight  gains  (5K)  Vs.y/ time  of  3502  epoxy 
coupons  subjected  to  various  stress  levels  under  constant  relative 
humidity  of  97%  and  at  temperature  of  40°C.  Predictions  of  Fick's 
law  are  drawn  by  dashed  lines.  After  Henson  and  Weitsman,  reference 
1371. 

Fig.  6:  Maximal  levels  of  moisture  weight  gain  {%)  attained  in  3502  epoxy 


samples  after  five  weeks  of  exposure  to  humidity  under  various  stress 
levels.  Conditions  same  as  in  Fig.  5.  After  Henson  and  Weitsman, 
reference  137). 
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Fig.  7;  Relative  moisture  absorption  (%)  in  epoxy  coupons  Vs.  time. 

Samples  immersed  in  distilled  water  at  5°C  under  various  levels  of 
compressive  stress.  After  Yaniv  and  Ishai,  reference  138]. 

Fig.  8:  Reduced  O2  permeability  in  extruded  PET  films  Vs.  stress.  After 
Seymour  and  Weinhold,  reference  1 39). 

Fig.  9:  Solubility  S  (in  %)  Vs.  non-dimensional  /time  of  epoxy  samples 
immersed  in  water.  Regime  "A"  -  sorption  behavior  at  75“C.  Regime 
"B"  -  increased  sorption  in  a  sample  equilibrated  at  75°C,  then 
cooled  to  23'C.  After  Apicella  et.  al.,  reference  136). 
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